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We present measurements of branching fractions and charge asymmetries for charmless B-meson 
decays to three-body final states of charged pions and kaons. The analysis uses 81.8 fb" 1 of data 
collected at the T(4S) resonance with the BABAR detector at the SLAC PEP-II asymmetric B 
Factory. We measure the branching fractions B(B + — > 7r + 7r~7r + ) = (10.9 ± 3.3 ± 1.6) x 10~ 6 , 
B(B + -> K+-K--K+) = (59.1±3.8±3.2)xl0~ 6 , and 23(B+ -> K + ' K~ K+) = (29.6±2.1±1.6) x 10~ 6 , 
and provide 90% C.L. upper limits for other decays. We observe no charge asymmetries for these 
modes. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The study of charmless hadronic B decays can make important contributions to the understanding of CP vi- 
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olation in the Standard Model, as well as to models 
of hadronic decays. Reference Q proposes that the 
interference between various charmless decays and the 
Xco resonance can be used to measure the Cabibbo- 
Kobayashi-Maskawa (CKM) angle 7, while the decay 
B + — > 7r + 7r~7r + can be used to reduce the uncertainties 
in the measurement of the CKM angle a 0] ■ We present 
branching fractions and charge asymmetries of charged- 
B-meson decays to three-body final states of charged pi- 
ons and kaons Q , with no assumptions about intermedi- 
ate resonances and with charm contributions subtracted, 
which allows us to set a tight bound on the charmless 
contribution to the measurement of 7 0. Upper limits 
and measurements of some of these branching fractions 
have been obtained previously with smaller statistics . 

The data used in this analysis were collected at the 
PEP-II asymmetric e + e~ storage ring with the BABAR 
detector, described in detail elsewhere The on- 

resonance data sample consists of 88.8 million BB pairs 
collected at the T(4S) resonance during the years 1999- 
2002. We also use 9.6 fb _1 of off-resonance data, col- 
lected 40 McV below the T(AS) resonance, to charac- 
terize the backgrounds from e + e~ annihilation into light 
qq pairs. We assume that the T(4S') decays equally to 
neutral and charged B meson pairs. 

Hadronic events are selected based on track multiplic- 
ity and event topology. Backgrounds from non-hadronic 
events are reduced by requiring the ratio of Fox- Wolfram 
moments H2/H0 [|| to be less than 0.98. Candidate B 
decays are formed by combining three charged tracks, 
where each track is required to have at least 12 hits in 
the tracking chamber, a minimum transverse momentum 
of 100 MeV/c, and to be consistent with originating from 
the beam-spot. 

Signal decays arc identified using two kinematic vari- 
ables: 1) the difference AE between the center-of-mass 
(CM) energy of the B candidate and \fs/2, where y/s 
is the total CM energy, and 2) the beam-energy substi- 
tuted mass toes = \J (s/2 + p t ■ Pb) 2 /Ef — p 2 B , where 
the B momentum p B and the four- momentum of the 
initial state (Ei,pi) are defined in the laboratory frame. 
The AE and toes distributions of signal events are Gaus- 
sian with resolutions of 20 MeV and 2.7 MeV/c 2 , respec- 
tively. The typical AE separation between modes that 
differ by substituting a kaon for a pion in the final state 
is 45 MeV, assuming the pion mass hypothesis. 

Charged pions and kaons are identified using energy 
loss (dE/dx) in the silicon detector and tracking cham- 
ber, and, for tracks with momenta above 700 MeV/ c, the 
Cherenkov angle and number of photons measured by the 
Cherenkov detector. The efficiency of selecting kaons is 
approximately 80%, which includes the geometrical ac- 
ceptance, while the probability of misidentifying pions 
as kaons is below 5%, up to a laboratory momentum 
of 4.0 GcV/c. Pions are required to fail both the kaon 
selection and an electron selection algorithm based on in- 



formation from dE/dx, shower shapes in the calorimeter 
and the ratio of the shower energy and track momentum. 

We remove B candidates when the invariant mass of 
the combination of any two of its daughter tracks (of 
opposite charge) is within 6a of the mass of the D° me- 
son or within 3a of the mass of the J/tpi ip(2S) or Xco 
mesons Q- Here, a is 10 MeV/c 2 for D°, 15 MeV/c 2 for 
J/tp and ^(25), and 18 MeV/c 2 for X c0- 

To suppress background from light-quark and charm 
continuum production, two event-shape variables are 
computed in the CM frame. The first is the cosine of 
the angle 9^ between the thrust axis of the selected B 
candidate and the thrust axis of the rest of the event, 
i.e. all charged tracks and neutral particles not assigned 
to the B candidate. For jet-like continuum backgrounds, 
|cos0^| is strongly peaked towards unity, while it is essen- 
tially uniform for signal events. For each signal mode we 
fix an upper limit on |cos0£|, between 0.575 and 0.850. 
This rejects between 95% and 85% of the background, 
depending on the decay mode. 

The second event-shape variable is a Fisher discrim- 
inant |||, which is formed from the summed scalar mo- 
menta of all charged and neutral particles from the rest of 
the event within nine 10°-wide nested cones coaxial with 
the thrust axis of the B candidate. The coefficients of 
the Fisher discriminant are chosen to maximize the sepa- 
ration between signal and continuum background events, 
and arc calculated for each signal mode separately us- 
ing Monte Carlo simulated signal and continuum events. 
A further 50% to 75% of the remaining background is 
rejected, depending on the decay mode, by applying se- 
lection requirements on this variable. 

B decay candidates passing the above selection criteria 
arc required to lie in a signal region defined as follows: 
|m ES - m B \ < 8 MeV/c 2 and \AE - (AE) | < 60 MeV, 
where (AE) = 7 MeV is the mean value of AE measured 
from on-resonance data for the calibration sample B~ — > 
D a n~,D° — > K~ir + , and tub is the nominal mass of 
the charged B meson Q. Figure Q shows the projections 
of the on-resonance data in the signal region onto the 
AE and toes axes. Each plot shows the expected levels 
of continuum and BB background, where the latter is 
parameterized from Monte Carlo samples. 

The residual continuum background level is estimated 
from the observed number of events in the grand sideband 
(GSB) region, defined to be 5.21 < m ES < 5.25 GcV/c 2 
and \AE - (AE) \ < 100 MeV, and extrapolating into 
the signal region. The shape of the toes distribution of 
the background is parameterized according to the phe- 
nomenologically motivated ARGUS function [|J, and is 
measured using off-resonance data and the upper side- 
band in the AE variable in on-resonance data (0.10 < 
AE < 0.25 GeV). A quadratic function is used to pa- 
rameterize the AE distribution of the background. The 
ratio of the integrals over the signal and GSB regions of 
the product of the AE and toes shape functions, R, gives 
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the ratio of the number of background events in the two 
areas. 
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FIG. 1: Projections of AE and mES for B + — » 7r + 7r 7t + 
(a and b), B+ -> ^+71^^+ (c and d), B+ -> A+A-7T+ 
(e and f), B+ -> A^A^A^ (g and h), B+ -> A~7r+7r + (i 
and j) and £> + — » A + A + 7r~ (k and 1) in the signal region. 
The signal region requirement was made on the orthogonal 
variable in each case. The dashed curves show the continuum 
background, while the solid lines include the BB background. 



The branching fraction for each channel is measured 
over the whole Dalitz plot, which is divided into 28 x 28 
cells of equal area (1 GeV 2 /c 4 ) 2 to enable us to find the 
selection efficiency as a function of position in the Dalitz 



plot. Taking €j to be the efficiency of reconstructing the 
signal in the i th bin of the Dalitz plot, determined from 
Monte Carlo simulated events, the branching fraction for 
each signal mode is given by: 



B 



N 



BB 



(N u - RN 2l - N xPi ) 



n x - rib , (1) 



where Nu and N 2 i are the number of events observed 
in the signal and GSB regions, respectively, while N x pi, 
n x and rib are background contributions that are defined 
below. N B g is the total number of BB pairs in the 
data sample. No significant differences were found for 
the value of R (defined earlier) in different regions of the 
Dalitz plot, so an average value is used for all bins. 

The probability of a kaon being misidentificd as a pion 
is 20%. This means there is significant cross-feed into 
the signal region from the decay mode that has one more 
kaon, which is subtracted for each bin, i. This is repre- 
sented by the N x pi/ei term in Eq. QJ, where N x is the 
total number of events that is the source of the cross-feed, 
and pi is the probability for the cross-feed events to pass 
the selection criteria for the i th bin, which is estimated 
from Monte Carlo samples. The B + — > K + K~K + mode 
has N x = 0, since it has no cross-feed backgrounds. For 
the other decays, N x is obtained by multiplying N BB 
by the branching fraction of the signal mode that has 
a kaon substituting a pion in the final state. There is 
also second-order cross-feed where either two kaons are 
misidentificd as pions (probability of 4%), or one of the 
pions is misidentified as a kaon (probability of 2%). This 
is represented by the n x term in Eq. 

Finally, the rib term represents the small number of 
other BB backgrounds that are subtracted. For all signal 
modes except B + — > K + K~K + , rib is obtained from the 
number of D° and D° candidates whose invariant mass is 
beyond the 6er range. For B + — > K + ir~ir + , there is also a 
contribution from B — > T)'(—> p°7)AT ± decays, which is 
estimated from the selection efficiency from Monte Carlo 
simulated decays, and the branching fraction quoted in 
Ref. 0. By using a mixture of Monte-Carlo-simulated 
charm and charmless decays, we found that there were 
no other significant BB backgrounds. 

We do not divide the Dalitz plot into cells for the 
Standard-Modcl-supprcssed modes B + — > K~-k + tt + and 
B + — > K + K + ir~ , and instead use the average values of 
the signal efficiency and cross-feed terms. 

The branching fraction results are summarized in Ta- 
ble |U where the first four rows show the total number 
of events in the signal and GSB regions, the average sig- 
nal efficiencies (e), and the values of R for each mode. 
The absolute efficiency variation across the Dalitz plot 
typically varies between ±2% and ±5% from (e). 

Rows A and B represent the total number of events 
and the amount of continuum background in the signal 
region, corrected for efficiency. The uncertainties for row 



6 



TABLE I: Branching fraction results for on-resonance data. The quantities and their uncertainties are explained in the text. 



Signal Mode 


± T ± 

7T 7T 7T 




A^A'+TT* 


A^A'+A* 


A' + 7T ± 7r ± 


A ± A ± 7T + 




1029 


1502 


733 


646 


494 


209 




5577 


5209 


4012 


1308 


3268 


1025 


(e)(%) 


12.7 ± 0.5 


12.8 ± 1.4 


13.9 ± 0.9 


14.9 ± 0.9 


18.5 ± 0.9 


15.3 ± 0.7 


R 


0.144 ± 0.003 


0.146 ±0.003 


0.150 ±0.003 


0.158 ±0.006 


0.155 ± 0.003 


0.157 ± 0.006 


A) Y.Nu/ei 


7597 ± 275 


11056 ±327 


5071 ±216 


4011 ± 182 


2670 ± 120 


1366 ± 94 


B) T RNii /e, 

J—" j £ ^ j J- 1.-1 1 j LJ 


5938 ± 94 ± 117 


5604 ±89 ±111 


4041 ± 72 ± 80 


1381 ± 46 ± 55 


2738 ± 48 ± 53 


1052 ± 33 ± 40 




474 ± 33 ± 40 


22 ± 1 ± 30 


671 ±15 ±59 






344 ± 31 


D) n x 




-189 ±34 


110 ± 128 




— 


53 ±5 


E) D° Bkgnd 


216 ±24 


268 ± 28 


47 ±6 




33 ±5 


31 ±5 


F) 77' A Bkgnd 




106 ± 30 










G) Signal Yield 


970 ± 291 ± 130 5246 ± 339 ± 127 202 ± 227 ± 163 2630 ± 188 ± 55 


-101 ± 129 ± 53 


-114 ± 100 ±51 




±22 ± 50 


±39 ± 247 


±16 ±9 


±12 ± 124 


±0±5 


±0±5 


B (xlO" B ) 


10.9 ± 3.3 ± 1.6 


59.1 ±3.8 ±3.2 


2.3 ± 2.6 ± 1.8 


29.6 ± 2.1 ± 1.6 


-1.1 ±1.5 ±0.6 


-1.3 ±1.1 ±0.6 


Significance (<r) 


5.7 


> 6 


1.1 


> 6 






90% C.L. 






< 6.3 




< 1.8 


< 1.3 



A come from the statistical errors in Nu, while the un- 
certainties for row B correspond to the statistical errors 
in Na, and the systematic errors from R, which arise 
from the limited statistics in the sideband region and off- 
resonance data. 

Row C shows the expected background from cross- 
feed events. The first and second uncertainties of these 
quantities represent the systematic errors in pi and N x , 
respectively, except for the channel B + — ► A' + A' + 7r~, 
where the uncertainty represents the average of the pi 
and N x contributions. The second-order cross-feed terms 
n x are shown in row D. Note that the n x term for 
B + — > K + ir~-K + is negative, which corrects for the 
B+ K+K-K+ cross-feed into B+ K + K~-k + , 
which in turn contributes to the cross-feed background 
for B + — * K + tt~7t + . 

Rows E and F show the expected backgrounds from D° 
and rj'K decays, which include the uncertainties from the 
selection efficiencies and the branching fractions of the 
background decays Q ■ The sum of these two rows gives 
the value of n& in Eq. 

Row G shows the signal yield, obtained by subtract- 
ing rows B to F from row A. The first uncertainty is 
the combination of the statistical errors of the number 
of events in the signal and GSB regions. The second 
uncertainty corresponds to the sum in quadrature of all 
the systematic errors from rows B to F. The third error is 
from the bin-by-bin uncertainty of the selection efficiency 
This is zero for B+ -> K-tt + tt + and B+ -> K+K+ir', 
where we only use the average efficiencies. The last un- 
certainty originates from global systematic errors in the 
signal efficiencies due to charged-particle tracking (0.8% 
per track), event-shape variable selections (1.0 to 2.5%) 
and from particle identification (1.4% and 1.0% per pion 
and kaon track, respectively). 

The next row in Table [I] shows the branching fraction 



results, where the first uncertainties are from the sta- 
tistical errors in the number of events, while the second 
uncertainties are the sum in quadrature of all systematic 
errors mentioned above. 

The significance of each branching fraction result, un- 
der the null hypothesis, is defined as the ratio of the 
signal yield to the total (statistical and systematic) un- 
certainty of the background in the signal region. We ob- 
serve significant signals for the modes B + — > 7r + 7r~7r + , 
B+ -»■ K+TT-TT+ and B+ -> K+ K~ K + ', and provide 
90% C.L. upper limits for the other channels, using the 
formalism in Ref . ^3 ■ The branching fraction of the con- 
trol sample B~ — > D°ir~,D° — » K~ir + , which has the 
same final state as B + — * A' + 7r~7r + , is measured to be 
(190 ± 3 ± 10) x 10~ 6 , which agrees with the average of 
published measurements (201 ± 20) x 10" 6 0- 

We have also measured the charge asymmetries for 
the modes with observed signals using a method sim- 
ilar to that used for the branching fraction measure- 
ments. The charge asymmetries are defined as A — 
(N- - N+)/(N~ + N+), where N~ (N+) is the signal 
yield for negatively (positively) charged B candidates, 
as defined by row G in Table HJ The normalisation fac- 
tor N B g cancels out in the asymmetry ratio, while the 
cross-feed and BB background contributions cancel in 
the asymmetry numerator. The measured charge asym- 
metries are A(B+ -> tt+tt-t: + ) = -0.39 ± 0.33 ± 0.12, 
A(B+ -> A:+7r-7r+) = 0.01 ± 0.07 ± 0.03 and A(B+ -> 
K+K-K+) = 0.02 ± 0.07 ± 0.03, where the first un- 
certainties are statistical and the second are systematic, 
which include the charge bias of the tracking and particle 
identification selection requirements (1%)- 

In summary, we have measured the branching fraction 
of B + — > 7r + 7r~7r + for the first time, which is smaller 
than that assumed in Ref. , and we have also observed 
the channels B+ -»■ AT + 7r-7r+ and B+ -> K + K~K+. 
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We observed no charge asymmetries in these decays. 
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